Original Article 

Environmental Health and Toxicology 
Volume: 27, Article ID: e2012008: 10 pages 
http://dx.doi.org/10.5620/eht.2012.27.e2012008 elSSN 2233-6567 



Open Access 



Modeling Human Exposure Levels to Airborne Volatile Organic 
Compounds by the Hebei Spirit Oil Spill 

Jong Ho Kim 1 , Byoung Kyu Kwak 2 f Mina Ha 3 f Hae-Kwan Cheong 4 , Jongheop Yi 2 

1 Supply Chain Logistics Consulting Team, Samsung SDS, Seoul; 2 School of Chemical and Biological Engineering, Institute of Chemical Process, College of 
Engineering, Seoul National University, Seoul; 3 Department of Preventive Medicine, Dankook University College of Medicine, Cheonan; 4 Department of Social 
and Preventive Medicine, Sung kyu nkwan University School of Medicine, Suwon, Korea 



Objectives: The goal was to model and quantify the atmospheric concentrations of volatile organic compounds (VOCs) as the result of the 
Hebei Spirit oil spill, and to predict whether the exposure levels were abnormally high or not. 

Methods: We developed a model for calculating the airborne concentration of VOCs that are produced in an oil spill accident. The model 
was applied to a practical situation, namely the Hebei Spirit oil spill. The accuracy of the model was verified by comparing the results with 
previous observation data. The concentrations were compared with the currently used air quality standards. 

Results: Evaporation was found to be 10- to 1,000-fold higher than the emissions produced from a surrounding industrial complex. The 
modeled concentrations for benzene failed to meet current labor environmental standards, and the concentration of benzene, toluene, ortho- 
meta- para-xylene were higher than the values specified by air quality standards and guideline values on the ocean. The concentrations of 
total VOCs were much higher than indoor environmental criteria for the entire Taean area for a few days. 
Conclusions: The extent of airborne exposure was clearly not the same as that for normal conditions. 
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INTRODUCTION 

Accidental oil spills are occurred frequently and can have 
severe adverse effects on the marine environment and 
human health. During the period 1970 to 2007, about 46 oil 
spills occurred per year and accounted for one hundred 
fifteen thousand tons of spilled oil. In addition, many oil 
spills typically occur near coastal areas [1] . 

The oil tanker Hebei Spirit collided with a crane barge on 
December 7, 2007, 8 km off the west coast of the Taean 
region in South Korea. The Hebei Spirit subsequently 
released an estimated 12,547 kL of crude oil. About 58,000 
citizens live in Taean, and more than 870,000 persons were 
involved in clean-up activities, including volunteers, army 
personnel, police, public servant, inhabitants, and specialized 
companies. It is likely that the clean-up workers and 
inhabitants, including elderly and weak people were exposed 
to high levels of evaporated volatile organic compounds 
(VOCs), and that this exposure likely had an impact on the 
health of people who were present at that time. In fact, 
adverse health effects as the result of exposure to VOCs 
have been reported in several studies of other oil spill 
accidents [2-5] . However, at the case of the Hebei Spirit oil 



spill accident, subsequent initial airborne exposure 
measurements were not carried out to confirm this, as has 
often been the case with other oil spill accidents. Without 
observation data regarding airborne VOCs concentrations, a 
modeling study is the only solution for estimating inhalation 
exposure which plays an important role in assessing possible 
adverse health effects. 

Several research groups had developed computational 
models for calculating the fate of spilled oil [6-18]. These 
models generally focused on the movement of an oil slick on 
the sea surface from the standpoint of weather-related 
processes such as evaporation, emulsification and 
dissolution. These models have proven to be useful for 
predicting the fate of an oil slick and changes in the 
physicochemical properties of the slick, and have affected 
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contingency plans and immediate protection strategy for 
protecting the coast line from contamination. However, our 
current understanding of the dispersion of evaporated oil 
components into the air remains incomplete, although air 
dispersion is an important factor in the inhalation exposure 
of toxic volatile oil components for humans. 

Herein, a modeling technique is proposed that permits the 
atmospheric concentrations of VOCs to be assessed after an 
oil spill. The modeling technique was developed by 
integrating the weathering model and the California Puff 
(CALPUFF) air dispersion model, and was applied to the 
Hebei Spirit oil spill. Our main aim was to assess the levels 
of VOCs that clean-up workers, as well as inhabitants, were 
exposed to and to predict the adverse effect of these levels 
on human health. This modeling study resulted in important 
information that promises to be useful in predicting adverse 
effects of such components on human health and in future 
guidelines for subsequent epidemiological studies. 



gravity- viscous formulation proposed previously by Drivas 
[6] and Hoult [7] and later modified by Mackay et al. [10]. 
The spreading rate is calculated as: 



dA 
dt 



4 

A 



(1) 



where A = the area of the slick (m 2 ), V = volume of the oil 
slick (m 3 ), Ki = a. constant with a default value of 150 s 1 and 
t- time (s). 

To solve equation 1 by a numerical method, an analytical 
solution was adopted for use as the initial conditions. Lehr et 
al. [11] developed an analytical formula that is suitable for 
estimating the initial size of a spill given the observed area of 
the spill. 
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MATEIALS AND METHODS 

The current section describes the available algorithms and 
data that are currently used to assess the evaporation rate of 
VOCs considering other weather-related processes, and 
illustrates the methods used to predict the concentration of 
evaporated compounds in air by the air dispersion model. The 
algorithms for evaporation with other weathering processes 
typically draw upon previously developed algorithms. The 
dispersion and concentration of VOCs in an atmospheric 
environment is predicted by employing a CALPUFF model, 
and various types of observation data (wind, current, aerial 
photographs, satellite images, etc) was introduced into the 
models. Each algorithm and method is described individually. 
The integration of the different algorithms and methods needs 
to reflect the fact that they interact with each other. 

1. Accident Scenario 

In the Hebei Spirit oil spill, three types of crude oil, United 
Arab Emirates Upper Zakum, Kuwait export crude and 
Iranian heavy crude, were released over a period of about 6 
hours. To calculate changes in the physicochemical 
properties of the spilled oil, it is necessary to define their 
initial properties. The averaged physicochemical properties 
for the three kinds of crude oil were used as the initial 
properties. The values for the average American Petroleum 
Institute gravity, Reid vapor pressure, density and the 
dynamic viscosity of initial crude oil properties were 31.4, 
46 kPa, 865 kg/m 3 and 56 cP, respectively. 



where po= density of the oil (g/m 3 ), p = density of the 
water (g/m 3 ) and W = the wind speed (m/s). The spreading 
process generally ceases at a terminal thickness of 0.1 mm 
for heavy crude oils [17]. The algorithm assumes a uniform 
thickness h (m), which can be calculated as: 



V_ 
A 



(3) 



111. Evaporation 



Evaporation is generally considered to be the primary 
process resulting in the loss of light ingredients during the 
first few hours of an oil spill [9] . To assess the evaporation 
rate, we applied the pseudo-component approach that has 
been described in previous studies [8]. The pseudo- 
component approach characterizes the mass of the individual 
hydrocarbon components evaporated, grouped by vapor 
pressure fraction. 



dm i 
dt 
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where mi = mass of the i th component (g), Pi - vapor 
pressure (atm), fi = fraction of spill which is constituent i, R 
= the gas constant (8.206E 5 atm m 3 / mol K), T = tempreature 
(K) and K2 = mass transfer coefficient (m/s) given by 
Mackay et al. [10] , expressed as: 



K 2 =0m9W- 0J8 D- 0Al S' 0 - 67 - 2.5X10" 3 W 0J8 (5) 



11. Oil Slick Spreading 

The rate of spreading of the oil slick is based on the 



where W - the wind speed (m/s), D = spill diameter (m) 
and Sc = the Schmidt number. 
To apply the pseudo-component approach for identifying 
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Table 1. Weight fraction and physicochemical properties of crude oil ingredients for the oil spill weathering model [8] 



No. 


Material 


Weight 
fraction 


Molecular Average 
weight density 
(g/moL) (kg/m 3 ) 


Average 
solubility 

(g/m 3 ) 




Vapor pressure (Pn, mm Hg) at a 
given temperature (T, K) 


1 


Paraffin (C6 - CI 2) 


0.15 


oO - 1 /U 


710 

/ 1U 


4.755 


Lo£ 


r p 1= 6.94-1417.61 /(-70.98+T) 


2 


Paraffin (CI 3 - C25) 


0.15 


184-352 


770 


0.007 


L0£ 


r P 7 = 7 01 -1825 04/(-123 39+T) 


3 


Cycloparaffin (C6 - CI 2) 


0.2 


84 - 164 


810 


28 


Lo£ 


I P3= 6.91 -1441. 79/(-68.45+T) 


4 


Cycloparaffin (CI 3 - C23) 


0.2 


156-318 


900 


0.5 


Lo£ 


I P4= 6.99-1893 .78/(-121.33+T) 


5 


Aromatic mono and dicyclic (C6 - CI 1) 


0.042 


78 - 143 


940 


890 


Lot 


I P 5 = 6.91 -1407.34/(-64.67+T) 


6 


Aromatic poly - cyclic (C 12 - CI 8) 


0.03 


128 - 234 


1000 


6.25 


Lo£ 


rP 6 =6.97-1801.00/(-110.38+T) 


7 


Naphtheno - aromatic (C9 - C25) 


0.07 


116-300 


980 


0.5 


L0£ 


I P 7 = 6.97-1789.85/(-108.59+T) 


8 


Residual including heterocycles 


0.15 


300 - 900 


1010 


0 


Ps= 


0 


9 


Benzene 


0.001 


78.11 


878.6 


800 


Lo£ 


I (P 9 /7.6)= -8.43Log(T)+(-6281/T)+71.11+6.2- 6 T 2 


10 


Toluene 


0.003 


92.14 


867 


526 


L0£ 


I (Pio/7.6)= -8.80Log(T)+(-6918/T)+74.14+5.75- 6 T 2 


11 


Ethylbenzene 


0.001 


106.17 


867 


206 


L0£ 


I (Pn/7.6)= -9.55Log(T)+-7638/T)+79.79+5.65- 6 T 2 


12 


ortho-, meta-, para-Xylene 


0.003 


106.17 


870 


0 


Pl2 = 


=5.16 



C, number of carbon atoms. 



the evaporation profile of chemicals, the physicochemical 
properties of the crude oil are required. The range of 
physical and chemical properties for each fraction is 
summarized in Table 1 . 

IV. Emulsification 

Under the influence of sea waves, water droplets are 
entrained by the oil slick, resulting in the formation of a 
water-in-oil emulsion. Mackay et al. [10] proposed the 
following expression for calculating the rate of incorporation 
of water into an oil slick: 



^ = 2.0X10" 6 (^ + 1) 2 (1-^) 
dt C l 



(6) 



where Yw = the fraction of water in the oil, Ci = final water 
content fraction (0.7 for crude oils and heavy fuel oil). 

V. Vertical Dispersion and Mechanical Removal 

Under the action of waves, especially breaking waves, the 
oil slick tends to break up, resulting in the formation of small 
oil droplets (diameter < 0.1 mm). The small droplets are 
incorporated in the water column during the natural 
horizontal dispersion process [18]. These processes are 
defined as natural vertical dispersion. 

In many oil spill accidents, dispersants and mechanical 
removal techniques are utilized to prevent long range 
horizontal dispersion and coastal pollution. The dispersant 
increases the velocity at which oil droplets are incorporated 
into the water column, and mechanical removal can prevent 
an oil slick from reaching a coastal area. In this study, the 
vertical dispersion (natural dispersion and accelerated 
dispersion by dispersants) and mechanical removal are 
regarded as total removal process, because these processes 
decrease the evaporation rate of volatile materials. The 
removal process can be expressed as follows: 



dV 
dt 



DV- DV- V m 



(J) 



where Dn and Dd = the entrainment coefficient cause of 
the natural dispersion and dispersant spraying, respectively 
(s -1 ), The Vm (m 3 /s) = the volume of oil components removed 
by mechanical removal. Reed [17] computed Dn by 
employing the entrainment formulation of Mackay et al. 
[10] , which is expressed as: 



D n -D a D h =[0.11(^ + l)2][l + 50/^ 2 MJ 



(8) 



where Da - the fraction of sea surface dispersed per hour, 
Db = the fraction of the dispersed oil not returning to the 
slick, W= wind speed (m/s), ju= viscosity (cp), St = oil- water 
interfacial tension (dyne/m). Dd and the Vm are obtained 
from the catalogues of the dispersants and published data in 
National Emergency Management Agency [19] . 

VI. Volume Change 

Volume change can be expressed as a function of the 
evaporation and removal process. In the case of volume 
change by an evaporation process, we assumed that the 
volume of the total mixture is the sum of the volume of the 
individual materials. Molecular interactions of individual 
ingredients in the mixture were ignored. 

dt tl dt Pi " d m (9) 

where Pi - density of i th component (kg/m 3 ). 

VII. Viscosity and Density 

An increase in viscosity due to mousse formation can be 
computed by the following equation [18]: 
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Figure 1. Geographical representation of the location of the oil spill and the area used in the simulation. 



dt 



c>2 

i-l 



2.5 g 

(1-C 3 7 r ) A 



(10) 



where C2= a constant between 1 and 10, where 1 denotes 
for light substances such as gasoline, and 10 for crude oils, Cj 
= final fraction water content (0.7 for crude oils and heavy 
fuel oil). The density increase may be calculated by [9]: 



= Y W P W + (1-Y W )(P 0 +C 3 F E ) 



VTI1. Transport on the Sea Surface 



(11) 



In order to calculate the horizontal dispersion of an oil 
slick on the sea surface, we analyzed satellite images and 
aerial photographs at specific times. The analyzed data were 
linked with the geographical information system (GIS). 
Finally the GIS -based snapshot of oil slick movement at a 
specific time was interpolated to estimate hourly motion. 

IX. Airborne VOC Concentration 

The CALPUFF model was integrated with the weathering 
algorithm to simulate the air dispersion of oil components. 
CALPUFF is one of the most widely used unsteady-state air 
dispersion models. This model is suitable for calculating the 
space and time of various air dispersed states. 

The location of the accident was approximately 8 km from 



the northwest of Malipo beach, the Taean district, Korea, 
235.6 km on the x coordinate and 4,084.51 km on the y 
coordinate of the 52 N universal transverse mercator (UTM) 
coordinate system (Figure 1). This area around the Taean 
district, which extends from 200 to 350 km on the x 
coordinate and 4,000 to 4,150 km on the y coordinate of the 
52 N UTM, was selected to run the model. To clearly 
discern the exposure level of volatile organic compounds in 
Taean, 57 numbered receptors were placed on the ocean 
floor, the coastline and the inland area. 

The Korean Meteorological Administration [20] provides 
free, hourly meteorological data such as wind speed, wind 
direction, temperature, rainfall, snowfall, cloud cover, solar 
radiation, and air pressure. In this study, we used the three 
meteorological data sets from three sites (described as a 
triangle in Figure 1); An-do, Ong-do, and Boryeong. Terrain 
map and landuse maps of the area were obtained from the 
environmental geographic information system [21]. The 
meteorological and geographical data were preprocessed 
(using the California Meteorological preprocess model 
(CALMET) meteorological preprocess model in the 
CALPUFF model) according to a literature report [22] and a 
previous study [23]. The CALMET model permits 
additional information to be calculated, such as the height of 
the Planetary Boundary Layer using the Monin-Obukhov 
similarity theory. Finally the CAPUFF model was run by 
preprocessed data. 
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Figure 2. Evaporation rate of spilled oil with time; (a) representative toxic volatile organic compounds; (b) light components of 
crude oil. 



RESULTS 

Figure 2 represents the total rates of evaporation of 
benzene, toluene, ethylbenzene and xylene (BTEX), and 
paraffin (C6-C12), cycloparaffin (C6-C12) and aromatic 
mono and dicyclic components (C6-C11). These results 
were obtained through applying the measured 
meteorological and hydrological situation for the Hebei 
Spirit oil spill to the model. The rates of evaporation 
fluctuate with time because evaporation is strongly 
influenced by the temperature and wind velocity. As 
depicted in Figure 2, the evaporation of benzene is almost 
complete in 10 hours, toluene, ethylbenzene and xylene have 
all nearly completely evaporated within 48 hours and the 
other volatile organic compounds had evaporated within a 
period of 4 days. The maximum emission rates of BTEX, 
and paraffin (C6-C12), cycloparaffin (C6-C12) and aromatic 
mono and dicyclic components (C6-C11) are 2,042, 3,115, 
697, 2,171 , 49,580, 62,549 and 16,612 in kg/hr, respectively. 

The tempo-spatial variation of atmospheric TVOCs 
[including BTEX and paraffin (C6-C12), cycloparaffin (C6- 
C12) and aromatic mono and dicyclic components (C6- 
Cll)] concentrations is illustrated in Figure 3. Until 7:00, 
December 10 (48 hours after the oil spill accident), the 
TVOCs moved in a southeasterly direction, following the 
coastline from the accident point. After 7:00, December 9, 
the downwind direction gradually changed to the west from 
the southeast. Finally, on 7:00 December 10 (78 hours after 
the oil spill), the evaporated components did not move 
inland. The maximum air concentrations of TVOCs are 
about 100,000 /g/m 3 on the ocean floor, 12,000 /g/m 3 at the 



coastline and the area of human settlement. 

Figures 4(a) - (c) illustrate the cumulative, maximum 1 hour 
averaged (1-hr) and maximum 8 hours averaged (8-hrs) 
concentrations for BTEX and TVOCs at each receptor point. 
The concentrations at 1-hour and 8-hours were selected from 
all of the recorded 1-hour and 8-hours averaged concentrations, 
therefore the data sampling times differed for chemical species 
and receptor sites. The receptor numbers can be broadly 
classified into 6 groups. Receptor number 0 represents the 
point of the oil spill on the ocean floor, receptors 1 to 4, 5 to 
1 1 , 12 to 22, 23 to 29 and 30 to 38 are located to the north, 
northwest, west, southwest and south side of the Taean 
district, respectively (Figures 1 and 3). The other receptors 
are situated in inland areas. 

On the ocean floor (receptor number 0), all cumulative 
VOCs concentrations are from 10-fold to 1 ,000-fold higher 
than the concentrations at the other sites. The cumulative 
concentrations of BTEX and TVOCs on this site were 
24,364 /g-h/m 3 , 42,716 /g-h/m 3 , 11,722 /g-h/m 3 , 33,978 /#g- 
h/m 3 and 589,571 /g-h/m 3 . The lowest concentrations of BTEX 
and TVOCs were modeled at receptor numbers 1, 2, 3 and 4, 
the north side of Taean. It is noteworthy that the respective 
concentrations for BTEX and TVOCs at receptor number 1 are 
about 10,000-fold, 2,000-fold, 150-fold, 800-fold and 100-fold 
lower than the concentrations at the coastline and inland areas. 
The cumulative concentrations of BTEX and TVOCs at 
receptor number 1 are 0.003 /g-h/m 3 , 0.187 /g-h/m 3 , 1.492 /g- 
h/m 3 , 0.511 /g-h/m 3 and 545.193 /g-h/m 3 . The concentrations 
of toluene, ethylbenzene, xylene and TVOCs at the coastline 
and inland are similar; the range of concentrations are 115.6 - 
574.4, 54 - 571.5, 117 - 783.6, and 9,585 - 162,463 /g-h/m 3 , 
respectively. In the case of benzene, the concentrations at the 
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Figure 3. Contour plot of the concentration of total volatile 
organic compounds by time and space. 

coastlines (2.2 - 106 /g-h/m 3 ) are, on average, 4 times higher 
than that at inland areas (1-28.3 /g-h/m 3 ). 

The spatial patterns for the concentration of toluene, 
ethylbenzene, xylene and TVOCs at 1 hour and 8 hours 
showed similar spatial patterns for cumulative 
concentrations. However, the cumulative concentrations for 
benzene are lower than ethylbenzene concentrations at each 
coast line, the 1-hour and 8-hours concentrations for benzene 
are similar to or higher than ethylbenzene concentrations at 
the coastline. This can be attributed to the high volatility of 
benzene and related meteorological conditions. 

DISCUSSION 

In early stage of an oil spill, several observations are 
usually carried out to confirm the airborne concentrations of 
BTEX. At 100 hours after the oil spill accident, the observed 
BTEX concentrations were 0.43-1.94 ppb, 0.58-1.59 ppb, 
0.16-0.65 ppb and 0.09-5.11 ppb, respectively [24]. The 
Ministry of Environment of Korea has been collecting data 
on these materials at monthly intervals. The annual average 
concentrations for these materials under normal conditions 
were 0.65 ppb, 1.158 ppb, 0.367 ppb and 0.499 ppb in the 
Taean area [21]. The difference in concentrations between 
normal conditions and the situation 100 hour after the oil 
spill were quite small. The reason for this is that the above 



four materials had completely evaporated in 48 hours, as 
shown in Figure 2(a), and because the wind direction 
changed to a westerly direction, the opposite direction from 
the Taean area by 78 hours after the oil spill (Figure 3). It 
necessarily follows that the observed concentration levels 
between normal conditions and 100 hours after the oil spill 
would be nearly the same. This result suggests that a national 
contingency plan for oil spill accidents must incorporate a 
monitoring plan that rapidly assesses the exposure of 
humans to VOCs. 

Initial observation data for the airborne concentration of 
volatile organic components is not available for the first 100 
hours, and the BTEX and TVOCs also evaporated very 
quickly. Therefore, it is not possible to directly confirm 
whether or not the modeled output from CALPUFF is 
accurate. This is the limitation of this study, which must 
contain precise values regarding concentration. So, we 
indirectly allowed some insight into the accuracy of the 
model with previously collected observation data taken from 
other oil spill cases. Jones et al. [25] performed a series of 
measurements on the concentration of airborne benzene over 
a water-filled moat at Little Sand Island in Mobile Bay in the 
US. In this study, almost all of the benzene had evaporated 
within about 8.5 hours. This result is consistent with our 
modeling results [Figure 2(a)] even if the environmental 
situations are different from the Hebei Spirit case. In 
addition, the benzene concentration above the oil slick was 
from a few thousand fjg/m 3 to few ten thousand fjg/m 3 for 8 
hours. This is very similar to our results for the ocean floor 
[Figure 4(a)] , in spite of the fact that the author's result tends 
to underestimate Jones's result. The underestimation is 
probably due to differences in meteorological and 
hydrological conditions, spill volume and receptor height 
(the receptor height of our modeling study was 1 .5 m and the 
observation height of the previous study was 19 mm above 
the slick). It can thus be concluded that the accuracy of the 
modeled evaporation rate of volatile components is quite 
acceptable. 

The emission rate for BTEX appears to be comparable 
with national data for emission rates [26,27]. The hourly 
emissions from an adjacent industrial complex (as depicted 
in Figure 1) for BTEX were, respectively, 7.7 kg/hr, 47.2 
kg/hr, 10.3 kg/hr and 38.1 kg/hr in 2007. The values for the 
maximum evaporation of BTEX are 265, 66, 3,319 and 56 
fold higher than the hourly industrial emissions of this area 
(Figure 2). This result clearly indicates that the airborne 
exposure level was abnormal, and supports the possibility 
that the clean-up workers and inhabitants living in adjacent 
areas were exposed to high levels of volatiles. 

The data provide solutions to the question of when and 
where attention should be focused in investigations of the 
impact of such situations on human health such as this one 
(Figure 3). Many of the clean-up workers and most of the 
inhabitants in Taean probably inhaled the large amounts of 
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Figure 4. Concentrations of benzene, toluene, ethylbenzene, xylene and total volatile organic compounds by the receptor 
number. 

a, cumulative concentrations; b, 1 hour averaged maximum concentrations; c, 8 hours averaged maximum concentrations. 



VOCs at the initial stage of the spill, because the VOCs were 
transported in a southeasterly direction from location of the 
accident for a period of 3 days. After 7:00, December 10 (72 
hours after the oil spill accident), the downwind direction 
changed to the west from the southeast direction. As a result, 
volatile components did not move toward the land. After this 
time, only clean-up workers on the ocean floor and at the 
coastal areas were exposed to excessive levels of VOCs. 
Clean-up workers on the ocean floor and on the coast line 
were exposed for more than 5 days, and inhabitants living in 
inland areas were exposed for 3 days. The population of the 
Taean area is about 64,000, and a total of 59,155 clean-up 
workers including inhabitants (9,887), police officers and 
army personnels (15,273), firemans (10,366), night guards 
(1,700), Civil Defense Corps members (2,663), the public 
officers (7,093), and volunteers (12,173) were involved in 
clean-up activities during the first 5 days. It would be 
expected that the exposure of humans to airborne volatile 
compounds might have an impact on the health of humans 
who live in the Taean coastal area, and especially clean-up 
workers who were exposed during the initial clean-up work 
(about 5 days). The inhabitants who live in the area north of 
Taean (receptor numbers 1,2,3 and 4 in Figure 1) were 



exposed to the lowest level of VOCs. These inhabitants 
could be used as a control group in investigating the health 
effects of volatile components. 

The previous database on environmental criteria provides 
information on risk. The concentration from the occurrence 
of the oil spill was very high compared with the normal 
conditions. So, we compared the concentration with 
workplace standards for each material, containing a high 
level standard. The Ministry of Labor of Korea [28] has 
regulated the extent of labor exposure by the workplace 
standards contained by the Korean Occupational Safety and 
Health Act (KOSHA). The KOSHA standards include 
threshold limit value - time weighted average (TLV-TWA) 
and threshold limit value - short term exposure limit (TLV- 
STEL). TLV-TWA is based on an exposure period of 8 
hours-per-day and a 40 hour-per- workweek schedule. TLV- 
STEL is the concentration that has been determined that 
workers can be continuously exposed for short periods of 
time. The TLV-TWA for benzene, toluene, ethylbenzene 
and xylene are 1 ppm, 50 ppm, 100 ppm and 100 ppm 
(approximately 3,000, 188,000, 435,000 and 435,000 in 
fjg/m 3 ), and the TLV-STEL are 5 ppm, 150 ppm, 125 ppm 
and 150 ppm (approximately 16,000, 560,000, 545,000, and 
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Table 2. Environmental air quality standards or guideline values with modeled maximum concentrations (various averaging 
time) 



Material 


Standards or guidelines (averaging time) 


Site 


Modeled maximum 
concentrations 

(jug/m 3 ) 


Averaging time 


Benzene 


3-5 fjg/m 3 (annual average) [31,37] 


Ocean 


2.78 


iy 






North 


3.5E-7-5.3E-5 








Coastline 


2.5E-4 - 0.012 








Inland 


1.1E-4- 0.003 






30 fjg/m 3 (1 h) [30] 


Ocean 


8,927.96 


lh 






North 


0.002 - 0.28 








Coastline 


1.21-52.3 








Inland 


0.48 - 13 




Toluene 


27-30 ^/m 3 (24 h) [31] 


Ocean 


1,774.99 


24 h 






North 


9.76-23.56 








Coastline 


0.94 - 22.49 








Inland 


0.007 - 24.26 




Ethylbenzene 


3800 fjg/m 3 (24 h) [33] 


Ocean 


466.51 


24 h 






North 


5.56-18.92 








Coastline 


0.49 - 14.53 








Inland 


0.006- 15.65 




Xylene 


250 ppb (1087.5 pg/m 3 ) (24 h) [32] 


Ocean 


2,799.35 


24 h 






North 


22.05 - 60.45 








Coastline 


1.97-53.95 




Total volatile 




Inland 


0.01-56.19 




organic 


300-500 pg/m 3 (2 h) [36] 


Ocean 


95,175.65 


2h 


compounds 




North 


137.91 -2,260.69 








Coastline 


1,387.92- 10,552.1 








Inland 


1,853.26-3,613.24 





655,000 in /g/m 3 ). As shown in Figure 2, the 1-hour maximum 
concentrations for BTEX were lower than the TLV-TWA. 
The 8-hour maximum concentrations for toluene, 
ethylbenzene and xylene were also lower than the TLV- 
STEL at all spatial positions. However, in the case of 
benzene on the ocean floor, the 1-hr maximum concentration 
(8,928 pg/m 3 ) is 56 % the level of TLV-STEL, and the 8- 
hours maximum concentration (3,037 fjg/m 3 ) is higher than 
that for TLV-TWA. The benzene concentration on the ocean 
floor failed to meet the workplace standards. These results 
corroborate the conclusion that the evaporated benzene 
could have had an adverse influence on the health of clean- 
up workers. Furthermore, the workplace standards were 
produced on the basis of normal healthy subjects; it has 
uncertainty in terms of assessing adverse effects on sensitive 
subgroups such as children, the elderly and pregnant women. 

Other available toxicity data are the ambient air quality 
standards or guidelines for VOCs. We illustrate the air 
quality standards or guidelines and modeled concentrations 
in Table 2. The Korean government recently established a 
standard value for benzene in ambient air of an annual- 
average of 5 //g/m 3 which corresponds with the European 
Union recommended annual limit. These values are around 
10 6 level in excess of carcinogenic risk based on the 
integrated risk information system [29] . The modeled annual 
average concentrations were 2.78 fjg/m 3 on the ocean and 
0.01 fjg/m 3 (maximum value) at site number 28. The annual 
average concentrations met the environmental standards and 



the guide line values (although the modeled values were 
similar to the Japanese standard). 

For the short-term criterion of air quality, Alberta 
Environment of Canada [30] published a guideline value for 
1-hr benzene concentrations (30 /g/m 3 ). As shown in Figure 
4(b), the maximum 1-hr benzene concentrations were 8,927 
fjg/m 3 on the ocean and 52.3, 49.95 and 49.95 fjg/m 3 at site 
numbers 20, 21 and 22, respectively. These values were 
much higher than the Alberta Environment guideline value. 
In the case of toluene, we compared our modeled data with 
the Swedish and WHO recommended air guidelines [31]. 
The modeled values at most of the receptor sites were lower 
than these guideline values; however the concentration on 
the ocean floor was 60 times higher than the guideline value. 
The Department of Environment and Heritage of the 
Austrian Government published guideline values for xylene 
[32]. The modeled xylene concentrations at coastal and 
inland areas were much lower than Austrian government 
guideline values. On the ocean floor, the modeled values 
were 2.57 fold higher than guideline values. 

The other available criterion is the guideline values for 
indoor air quality. These guideline values were extensively 
arranged [33]. The guideline value for ethylbenzene is 3,800 
fjg/m 3 , therefore the modeled ethylbenzene concentrations at 
all locations in Taean met the indoor air quality guideline 
value. The concentrations of TVOCs can be compared with 
the indoor air quality guideline values of TVOCs, although the 
ingredients and their portions are probably different from our 
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TVOCs. As shown in Table 2, most of the modeled 
concentrations did not meet the guideline values at most of the 
sites in Taean, except for receptor number 1 (137.91 /g/m 3 ). In 
a comparison study between the modeled results and air 
quality standards or guideline values, it can be concluded that 
the modeled concentration was clearly not within the 
acceptable guidelines of environmental quality. This suggests 
that the risk levels are probably quite high and it is possible 
that this had an adverse impact on human health. Therefore, 
follow-up studies were needed for the evaluation of the risk by 
BTEX and TVOCs. 

The published short term epidemiological studies are 
consistent with this conclusion. Using our model outputs, 
several epidemiological studies were carried out, to evaluate 
the acute health effect on pregnant women and children, and 
the results have been published [24,34]. In these studies, 
pregnant women exposed to higher ambient cumulative 
levels of BTEX were significantly more likely to report skin- 
related symptoms in the first day after the accident, and 
significantly more likely to report abdominal pain during the 
first through fourth days following the accident. Children 
who were exposed to high levels of BTEX reported more 
physical and mental symptoms than those who lived farther 
from the polluted coast line. These studies suggest that the 
exposure level to BTEX was abnormally high and is 
correlated with an increased risk of health effects among 
pregnant women and children. These results indirectly show 
that the modeled output has an acceptable accuracy for the 
tempo-spatial patterns of VOCs concentrations. 

In this study, an integrated procedure for modeling oil 
slick weathering and atmospheric dispersion is proposed, in 
an attempt to understand the exposure level of volatile 
compounds in a crude oil spill. We conclude that our 
modeling study has the potential to enhance current 
exposure analysis techniques, that they are feasible, and can 
provide baseline information for future long term 
epidemiological studies. In the case of an oil spill near a 
heavily populated area (such as the Hebei Spirit oil spill), a 
modeling study of airborne exposure is significant. 

From the labor environmental standards and air quality 
criteria, the findings showed that the exposure levels to 
clean-up workers and inhabitants were quite high. The 
results show that the adverse impact on humans is not only a 
possible short term effect but also may well persist for a 
longer term. Furthermore, future possible long-term impacts 
are unpredictable, because the available scientific data 
cannot predict the all possible adverse effects. Crude oil 
contains many types of hazardous compounds; therefore it is 
not possible to predict possible synergic effects between 
different compounds due to our current lack of knowledge. 
Therefore long term studies will be required to assess the 
adverse effect of events such as these on human health in the 
Taean area. 

Despite the advantages of the modeling study, it is 



important to seek further improvements. Future studies will 
be required to verify the proposed modeling method, 
including a comparison study of model outputs with 
observation data on oil spill situations and an uncertainty 
analysis. Additional future work should be undertaken to 
simulate the long term exposure of concentrations of weak 
volatile compounds (e.g. polycyclic aromatic hydrocarbons, 
persistent organic pollutants) in the air. 
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